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We study the finite temperature behavior of light scalar and pseudoscalar meson proper- 
ties in the context of a three-flavor nonlocal chiral quark model. The model includes mixing 
with active strangeness degrees of freedom, and takes care of the effect of gauge interactions 
by coupling the quarks with the Polyakov loop. We analyze the chiral restoration and decon- 
finement transitions, as well as the temperature dependence of meson masses, mixing angles 
and decay constants. The critical temperature is found to be T c ~ 202 MeV, in better agree- 
ment with lattice results than the value recently obtained in the local SU(3) PNJL model. 
It is seen that above T c pseudoscalar meson masses get increased, becoming degenerate with 
the masses of their chiral partners. The temperatures at which this matching occurs depend 
on the strange quark composition of the corresponding mesons. The topological susceptibil- 
ity shows a sharp decrease after the chiral transition, signalling the vanishing of the U(f )^ 
anomaly for large temperatures. 
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I. INTRODUCTION 



The detailed understanding of the behavior of strongly interacting matter under extreme con- 
ditions of temperature and/or density has become an issue of great interest in recent years. In 
this context, it is clearly important to study how hadron properties (masses, mixing angles, decay 
constants, etc.) get modified when hadrons propagate in a hot and/or dense medium. In partic- 
ular, since the origin of the light scalar and pseudoscalar mesons is related to the phenomenon of 
chiral symmetry breaking, the temperature and/or density behavior of their properties is expected 
to provide relevant information about a possible chiral symmetry restoration. Unfortunately, even 
if a significant progress has been made on the development of ab initio calculations such as lattice 
these are not yet able to provide a full understanding of the QCD phase diagram and 
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the related hadron properties, due to the well-known difficulties of dealing with small current quark 
masses and finite chemical potentials. Thus it is important to develop effective models that show 
consistency with lattice results and can be extrapolated into regions not accessible by lattice cal- 
culation techniques. In previous works [J-^Lthe study of the phase diagram of SU(2) chiral quark 
models that include nonlocal interactions [8j has been undertaken. These theories can be viewed 
as nonlocal extensions of the widely studied Nambu— Jona-Lasinio model [9]. In fact, nonlocality 
arises naturally in the context of several successful approaches to low-energy quark dynamics as, for 
example, the instanton liquid model [10] and the Schwinger-Dyson resummation techniques fill ]. 
Lattice QCD calculations 113] also indicate that quark interactions should act over a certain range 

n n 

in momentum space. Moreover, several studies [13|, [lj] have shown that nonlocal chiral quark 
models provide a satisfactory description of hadron properties at zero temperature and density. 
On the other hand, when looking at the description of the chiral phase transition, it has been 
noticed that for zero chemical potential these models lead to a rather low critical temperature T c 
in comparison with lattice results P, [5]. However, it has been recently shown that the inclusion 
of the Polyakov loop, which can be taken as an order parameter for the deconfinement transition, 
leads to a significant increase of the chiral restoration temperature both in two-flavor [15] and 
three-flavor [161 ] nonlocal models. The inclusion of the Polyakov loo p h as also been considered in 



the context of NJL-like models, namely the so-called PNJL models 
model [22]. 
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211 ] . and the quark-meson 



The aim of the present work is to go one step beyond previous analyses, studying the finite 
temperature behavior of light scalar and pseudoscalar meson properties in the context of three- 
flavor nonlocal chiral models that include mixing with active strangeness degrees of freedom, and 
taking care of the effect of gauge interactions by coupling the quarks with a background color gauge 
field. 

This article in organized as follows. In Sect. II we present the general formalism and derive the 
expressions needed to evaluate the different meson properties at finite temperature. In Sect. Ill 
we provide details concerning the determination of model parameter values as well as the results 
obtained at zero temperature. Our results for the behavior of the different meson properties as a 
function of the temperature are presented and discussed in Sect. IV. Finally, in Sect. V we sketch 
our conclusions. 
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II. THE FORMALISM 

We deal here with a nonlocal covariant SU(3) quark model which includes the coupling to a 
background color gauge field. The Euclidean effective action for the quark sector of this model is 
given by 

Se = J d 4 x U(x) [-ij„ D n + ™] 1>W ~ § [if (*) if (*) + if (*) J« (*)] 

-| T a6c [if (z)jg(x)jf!(x) - 3 if (*)] + W [A(x)]| , (1) 

where the chiral U(3) vector tp includes the light quark fields, ip = (u d s) T , and m = 
&\&g{m u ,md,m s ) stands for the current quark mass matrix. For simplicity we consider the isospin 
symmetry limit, in which m u = = rh. The fermion kinetic term includes a covariant derivative 
= d^ — iA^, where are color gauge fields, and the operator "J^d^ in Euclidean space is defined 
as X)i=i 3~fiW~ + 74^' w ith 74 = ^7o- Regarding the interaction terms, the currents ja' P (x) are 
given by 

if (x) = J d 4 z g(z) i, (x + |) A a tf, (x - ~) (2) 
£{x) = Jd*zg(z) + i-Y5X a ip(x-~) , (3) 

where g(z) is a form factor responsible for the nonlocal character of the interaction, and the matrices 
A a , with a = 0, .., 8, are the standard eight Gell-Mann 3x3 matrices — generators of SU(3) — plus 
Ao = \/2/3 13x3- The constants T a b c in the t'Hooft term (responsible for flavor-mixing) are defined 
by 

Tabc = ^ Cijk Cmnl ( A a ) im (Aft) jn (X c ) kl ■ (4) 

Finally, the action ([1]) also includes an effective potential IA that accounts for gauge field self- 
interactions. 

The partition function associated with the effective action Eq. ([T]) can be bosonized in the usual 
way introducing the scalar and pseudoscalar meson fields a a (x) and TT a (x) respectively, together 
with auxiliary fields S a (x) and P a {x). To deal with these auxiliary fields we follow the standard 
stationary phase approximation, which provides a set of equations that relate them to the scalar 
and pseudoscalar meson fields. Since we are interested in studying the behavior of various meson 
properties in the presence of a heat bath, we have to extend the bosonized effective action to finite 
temperature. In the present work this is done by using the Matsubara formalism. 
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The coupling of fermions to the color gauge fields is implemented through the covariant deriva- 
tive in the fermion kinetic term As usual, we will assume that the quarks move in a 
constant background field = zAq = ig 5^q GaX a /2, where Ga are the SU(3) color gauge fields. 
Then the traced Polyakov loop, which is taken as order parameter of confinement, is given by 
$ = |Tr exp(i/3</>), where (3 = 1/T, <fi = iAq. We will work in the so-called Polyakov gauge, in 
which the matrix <j> is given a diagonal representation <j) = (P3X3 + (ps^s, which leaves only two 
independent variables, ^3 and <f>%. 

To treat the resulting finite temperature system of interacting mesons in the presence of the 
Polyakov loop we consider first the mean field approximation (MFA), keeping only the nonzero 
vacuum expectation values a a . Note that due to charge conservation only cx a =o,3,8 can be different 
from zero. Moreover, 03 also vanishes in the isospin limit. The corresponding MFA grand canonical 
thermodynamical potential reads 



^MFA(T) = - 2 J2 [ Tr ln Vnc + £/(Pnc)] 
f,c Jp ' n 

+ U(*,T) , (5) 



f 



where / = u,d, s, c = r,g,b, and the shorthand notation J pn = Yin I d 3 p/(27t) 3 has been used. 
We have also introduced the definition p nc = {p ,u n — C ), where w n stands for the fermionic 
Matsubara frequencies and the quantities <fi c are defined by the relation <p = diag(0 r , (f) g , </!>&). The 
quark constituent masses T,f(p nc ) are here momentum-dependent quantities, given by 

s /(Pnc) = mf + a f g(p nc ) , (6) 

where g{p) is the Fourier transform of the form factor g(z). For convenience we have introduced 
mean field values cry given by 

/2 l [2 1 [2 2 

-00 + 03 + ^=^8 , a d = ^-a -a 3 + -j=a 8 , ° s = V 3 ^° + ^/f^ 8 ' 

and similar definitions hold for Sf in terms of So, S3 and S%. Note that in the isospin limit a u = a<i, 
thus we have T, u (p nc ) = S^(p nc ). Within the stationary phase approximation, the mean field values 
of the auxiliary fields Sf turn out to be related with the mean field values of the scalar fields df 



by 



23] 



a u + GS u + jS u S s = , a s + GS s + ^S 2 u = . (7) 
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The effective potential U(&,T), which accounts for Polyakov loop dynamics, can be fitted by 
taking into account group theory constraints together with lattice results, from which one can 
estimate the temperature dependence. Following Ref. 
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U{$>,T) 



211 ] we take 
a(T) $ 2 + b(T) m(l - 6 $ 2 + 8 <£» 3 - 3 $ 4 ) 



(8) 



with the corresponding definitions of a(T) and b(T). Owing to the charge conjugation properties 
of the QCD Lagrangian [24], the mean field traced Polyakov loop field <I> is expected to be a real 
quantity. Assuming that 4>z and 4>g are real- valued |2l|], this implies (f>8 = 0, = [2 cos(^3/T) + l]/3. 

For finite current quark masses the quark contribution to Qmfa(T) turns out to be divergent. 
To regularize it we follow the same prescription as in previous works [(J. Namely, we subtract 
from Qmfa(T) the quark contribution in the absence of fermion interactions, and then we add it 
in a regularized form, i.e. after the subtraction of an infinite, T-independent contribution. From 
the minimization of this regularized thermodynamical potential it is possible now to obtain a set 
of three coupled "gap" equations that determine the mean field values a u , a s and 4>s at a given 
temperature: 

<90mfa 



. 



(9) 



(<9o- u ,dcr s ,<9</> 3 ) 

In order to obtain the meson mass spectrum and other properties one has to consider the 
mesonic fluctuations around the mean field values. We begin by introducing a more convenient 
basis defined by 



where £ a = a a ,ir a , while i,j run from 1 to 3 (neutral fields are shifted by £ a 
scalar fields one has in this way 



(10) 

£„ - 0. For the 



/ ap (Js <?o 
V2 V6 V3 



\ 



\ 



V2 y/6 V3 

_o 2og (Jo_ 



(11) 



I J 



while a similar expression holds for the pseudoscalar sector, replacing ao — > vr, a — > rj and k — > 
K. Using this notation the resulting quadratic contribution to the finite temperature bosonized 
effective action can be written as 



S 



quad 



q,m 



(12) 
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where q m = (q ,v m ), v m = 2mirT being bosonic Matsubara frequencies. The functions Gf- kl in 
Eq. (fT2"j) are given by 

Gt kl (q 2 ,v 2 J = Ct 3 {q 2 y m ) Mi* + ((r*)- 1 )^ , (13) 

where 

\2 Pnc Pne ' ~F ^i(Pnc + Qm)^ 1 j {Pnc) 



Di(Pnc + qm)Dj(p 

r, 



(14) 



and 

rr 

r ij,fc« = G S H S jk ± y e ikh Cjlh S h ■ (15) 

In Eq. (Q3D we have defined Dj(s) = s 2 + S|(s), where j = 1, 2, 3 correspond to / = u, d, s in the 
notation of Eq. ©. 

Prom the quadratic effective action 5^ uad it is possible to derive the scalar and pseudoscalar 
meson masses as well as the quark-meson couplings. In what follows we will consider explicitly 
only the case of pseudoscalar mesons. The corresponding expressions for the scalar sector are 
completely equivalent, just replacing upper indices "— " by "+". In terms of the physical fields, 
the contribution of pseudoscalar mesons to S'g Uad can be written as 



oquad 
E 



p = \i {G n (q 2 ,vi) [n (q m ) 7T°(-<?m) + 2 n + (q m ) ^{-q^} 



q,m 



+G K {q 2 ) v 2 m ) [2K°(q m ) K°(-q m ) + 2 K + (q m ) K-(-q m j\ 
+G r) {q 2 ,vl l ) rj(q m ) 7](-q m ) + G^{q 2 ,u^) rf{q m ) rf{-q m )j ■ (16) 

Here, the fields rj and rf are related to the U(3) states 770 and % according to 

rj = cos 9 V rjs — sin Q n rjo 

r/' = sin^/ r/g + cos^/ % , (17) 

where the mixing angles are defined in such a way that there is no rj — rj mixing at the level 
of the quadratic action. In a similar way, in the scalar sector one has two physical scalar mesons 
a and /o(980) that are linear combinations of the states <7g and <tq, with mixing angles 6 a and 9f . 
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The functions Gp(q 2 , i^J introduced in Eq. (fT6j) are given by 



G n (q 2 , 



where 



and 



(G+-5 s )- i + C- u (g^^ 



(G+ -5 u )- i + C-(g^^) 



(18) 



G 88 (q 2 ,vL) + G 00 (q 2 ,V 



7t1 ,,2 1 



[Go 8 (g 2 ,^)] 2 + 



1 2 



G 88 {q 2 ,i / m) — G f oo(^' 2 > ^m) 



G 88 (q 2 , Vm) 
G 08 (q 2 , v^a) 



1 6G -4HS U -2HS S 
3 2G 2 -GHS S - H 2 S 2 
y/2 H{S S - S u ) 
3 2G 2 -GHS S - H 2 Sl 
1 6G + 4HS U - HS S 



+ C 88 {q 2 ,v 2 m ) 



+ C 08 (q 2 ,vl l 



3 2G 2 - GHSe - H 2 S 2 



Gq 8 — g {p uu C ss ) 



+ C 00 (q 2 y m 



a 



2 C«n + G ss 



00 



Now the pseudoscalar meson masses are obtained by solving the equations 



(19) 



(20) 



G P (-m 2 p,0) = 



(21) 



with P = it, K, r] and rf. The mass values determined by these equations correspond to the 
spatial "screening-masses" of the mesons' zeroth Matsubara modes, and their inverses describe 
the persistence lengths of these modes at equilibrium with the heat bath. It is worth to notice 
that there is a screening mass for each Matsubara mode. The full bound state propagator can be 
calculated via any polarization tensor that receives a contribution from the bound state, but only 
once all screening masses have been determined. The propagator obtained in this way is defined 
only on a discrete set of points along what might be called the imaginary-energy axis, and the 
"pole-mass", i.e., the mass that yields the bound state energy pole for q ~ 0, is obtained only 
after an analytic continuation of the propagator onto the real-energy axis. The fact that Lorentz 
invar iance is broken for T_> means that, in general, the pole mass and screening masses are 
not equal (see e.g. Ref. 



25]). Although the analytic continuation involved in this process is not 
unique, an unambiguous result is obtained by requiring that it yields a function that is bounded 
at complex-infinity and analytic off the real axis From this description it is nonetheless clear 
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that the screening masses completely specify the properties of T > bound states. The masses 
associated to the zeroth Matsubara mode studied here are spatial screening masses corresponding 
to a behavior exp(— mpr) in the conjugate 3-space coordinate r, and should correspond to the 
lowest state in each meson channel. In fact, these are the quantities usually studied in lattice 
calculations [3]. 

In the f] — rj' system, once the meson masses have been determined one can find the mixing 
angles 8^ and 9„/, which are in general different from each other. These are given by 

tan 29 P = ^""tf 2 m > P = V, v' ■ (22) 

G 00 (-m 2 p ,0) - G 88 {-m 2 p ,0) 

The meson fields have to be renormalized, so that the residues of the corresponding propagators 
at the meson poles are set equal to one. The corresponding wave function renormalization constants 
Zp are given by 



Z7} 



dG P (q 2 ,0) 



p dq 2 



(23) 



with P = tt, K, rj and rf . Finally, the meson-quark coupling constants Gp q are given by the 
original residues of the meson propagators at the corresponding poles, 

G Pq = Zp 12 . (24) 

In the case of the pseudoscalar mesons other important features are the corresponding weak 
decay constants / a &, defined by 

(0\Al(0)\7r b (q))=if ab q„. (25) 

where is the a-component of the axial current. For a, b = 1 . . . 7, the constants f ab can be 
written as 5 a b fp, with P = tt for a = 1, 2, 3 and P = K for a = 4 to 7. In contrast, as occurs with 
the meson masses, the decay constants become mixed in the a = 0, 8 sector. Details on how to 
obtain the expressions for the axial currents in the presence of nonlocal fields can be found e.g. in 
Refs. I3I, 23]. After a rather lengthy calculation we find that, at finite temperature, the pion and 



kaon decay constants are given by 

fn = 4/ uu (-m2,0) Z l J 2 , 



Sk = 2 [f us (-m 2 K , 0) + fsu(-m 2 K ,0)] zf , (26) 
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where 

Pnc " Qm^i iPnc + Qm) ~ (Pnc ' Qm + Pnc)^ j {Pnc) 



fij(q 2 ,v m ) =^2 \g(pnc + 

c J P> n I 



2 / Di (p nc + q m )Dj (p nc ) 

^i{Pnc + 



\2g [p nc +'y)-9 (Pnc) ~ 9 One + 



fa' 

2 - 



+ 9 [ Pnc + ~y 



Si [Pnc + ^pj + £j (pnc + ^) - St {Pnc) ~ (p nc + g r , 



(P 

Di{p nc + q m )Dj(p 

ncj 

In the case of the r] — r( system, two decay constants can be defined for each component (a = or 
8) of the axial current. They can be written in terms of the decay constants f a b and the previously 
defined mixing angles 9 V:V ' as 

= [/«8(-mJ J 0)cose„-/ a0 (-m;,0)sine„] Z}/ 2 
f% = [/«8(-mJ,,0)Hin^ + / a0 (-mJ,0)cos^] Z X J 2 . (28) 

Within our model, the decay constants f a b for a, b = 0, 8 are related to the defined in Eq. (|27p 
by 

/88(g 2 , ^m) = ^ [2/ ss (g 2 , Z^n) + fuu(q 2 , V m )] 

/oo(<f 2 ,^m) = g [2/ u «(g 2 ,f m ) + /ss(^ 2 ,^m)] (29) 

4\/2 

fo8(q 2 ,v m ) = fw{q 2 ,v m ) = — — [fuu{q 2 ,v m ) - f ss (q 2 ,V m )] . 

As expected, both the nondiagonal decay constants /os, /so and the mixing angles 9 n , 9^ vanish 
in the SU(3) symmetry limit. 

III. MODEL PARAMETERS AND ZERO TEMPERATURE RESULTS 

In this section we determine the model parameters to be used in our numerical calculations, and 
quote the results obtained for various meson properties at zero temperature. The latter include the 
values of meson masses, decay constants and mixing angles, as well as quark constituent masses, 
quark condensates and quark-meson couplings. 

At low temperatures, the value of the traced Polyakov loop is essentially determined by the 
effective potential in Eq. 0, therefore for T — > one has $ — > 0, cos(4>s/T) — > —1/2. Since for low 
T the Matsubara sums in the thermodynamical potential are governed by modes with large n, one 
has oj n — 4> c = [(2n + 1)tt — cfr c /T]T ~ uj n , thus for T — > the coupling of fermions to the Polyakov 
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loop vanishes. In this way, the zero-T calculations are similar to those carried out in Ref. |23|], in 
which SU(3) nonlocal chiral quark models without the inclusion of the Polyakov loop have been 
considered. As in that work, our numerical analysis has been performed using a Gaussian form 
factor, namely 

g(p) = exp ( V/A 2 ) , (30) 

which has been often considered in the literature. This introduces a new free parameter A, which 
plays the role of an ultraviolet cut-off momentum scale (we recall that the form factor is defined 
in Euclidean momentum space). At T = 0, the main difference between our analysis and that 
in Ref. [3] is that here we are considering an OGE-motivated nonlocal interaction, whereas in 
the previous work a different form [motivated by instanton liquid models (ILM)] for the nonlocal 
currents has been chosen. In the case of two-flavor models, a detailed comparison between these 



281 ]. showing that the results for both 



different interaction forms has been carried out in Ref. 

pi 

models are qualitatively similar. Notice that in Ref. [23j only the pseudoscalar meson sector was 
addressed. 

After the assumption of the form factor in Eq. (|30p . the nonlocal chiral quark model under 
consideration includes five free parameters, namely the current quark masses fh and m s , the cou- 
pling constants G and H and the cut-off scale A. In our numerical calculations we have chosen to 
fix the value of fh, whereas the remaining four parameters are determined by requiring that the 
model reproduces correctly the measured values of four physical quantities at zero temperature. 
These are the masses of the pion, kaon and rj' pseudoscalar mesons, and the pion decay constant 
/„-. Taking fh = 5 MeV, we obtain the following set of parameters: 

fh = 5 MeV (input) 

m s = 119 MeV 

A = 843 MeV 

GA 2 = 13.35 

HA 5 = -273.7 (31) 

Our numerical results are presented in Table I. For comparison, in the last column of this table 
we quote the measured values of meson masses, and the ranges in which the decay constants and 
mixing angles should fall according to most popular phenomenological approaches. Entries marked 
with an asterisk are those that we have taken as input values to fix the model parameters. From 
Table I it is seen that in general there is a reasonable agreement between the predicted meson 
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29l | . In addition, the 



masses and the empirical values quoted by the Review of Particle Physics 
obtained mass ratio m s /m = 23.8 is close to the corresponding current algebra prediction, namely 
m s /m = (2m 2 K — m^)/m 2 ~ 25. We point out that in the case of the mass of the k scalar 
meson the equation G K {— x 2 ,0) = has no solution in the real x axis. Hence we have defined 
the mass m K as the point where the absolute value of G R {— x 2 ,0) becomes minimal. A more 
sophisticated definition could be done by extending x to the complex plane, thus introducing a 
finite k width. In any case, this would not change significantly the mass value. A detailed analysis 
of the regularization prescriptions for the evaluation of loop integrals like those in Eqs. f)14[) and 
(|27l) has been carried out in Ref. [231 ]. 

Concerning the pseudoscalar meson decay constants, we notice that the predicted value for 
Jk is also phenomenologically acceptable. In fact, it turns out to be significantly better than 
that obtained in the standard NJL model, where the kaon and pion decay constants are found to 
be approximately equal to each other [9] in contrast with experimental evidence. Regarding the 
mixing angles and decay constants for the 77s — 770 system, the problem of defining and (indirectly) 
fitting these parameters has been revisited severe, thnes in the iiteratnre (see e„ Ref. Q, and 
references therein). As stated in the previous section, in general one has to deal with two different 
state mixing angles » P and foar decay constants ff,, where P _ „, ,< and « = 0, 8. This means that 
77 and 77' states do not need to be orthogonal, and the same occurs with /®,) and (/°, /°,) [3ll — 



33] . For the sake of comparison with phenomenological values of these parameters, we follow here 



Ref. [3l|] and express the four decays constants fp in terms of two decay constants f a and two 
mixing angles 9 a , where a = 0,8: 



fr, fr,\_lh cos 9 8 -fo sin6> 
frf /?7 \f8sm6 8 fo cos 9 



(32) 



In our framework the decay constants fp can be calculated from Eqs. (I28p . As shown in Table I, 
both the values obtained for #o as well as those obtained for the decay constants in the 77 — 77' 
sector are in agreement with phenomenological results. These have been taken from the analysis in 
Ref. [30] , in which the values obtained from different parameterizations have been translated to the 
four-parameter decay constant scheme given by Eq. (I32p . Notice that #8 and 6q are significantly 
different to each other, as occurs with the mixing angles 9„ and [which can be calculated from 
Eq. (I22h ]. This is in agreement with the analysis in Ref. 3l|], carried out within next-to-leading 



order Chiral Perturbation Theory and large Nc, which leads to 6% = —20.5°, 6q = —4°. 
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[ MeV ; 


5* 


(3.4 - 7.4) 


m s 


[ MeV ; 


119 


(108 - 209) 


Win 


[ MeV ; 


139* 


139 


m K 


[ MeV ; 


495* 


495 


m n 


[ MeV ; 


523 


547 


m v > 


[ MeV ; 


958* 


958 


m ao 


[ MeV ; 


900 


980 


m K 


[ MeV ; 


1380 


1425 


m„ 


[ MeV ; 


566 


400-1200 


m fo 


[ MeV ; 


1280 


980 


G-jzq 




3.98 




GKq 




4.30 




G m 




3.93 




G n 'q 




2.83 




@n 




-2.3° 




Or,' 




-40.3° 




o s 




-24° 


-(22° - 19°) 


Oo 




-7.7° 


-(10° - 0°) 


u 


[MeV; 


92.4* 


92.4 






1.17 


1.22 


fr\l fir 




1.14 


(1.17-1.22) 


frit fir 




0.16 


(0.11-0.37) 


frj' 1 fir 




-0.49 


-(0.42-0.46) 


frj'l /tt 




1.16 


(0.98-1.16) 



(*) Input values 

TABLE I: T = model predictions for various meson properties: masses, mixing angles, decay constants 
and quark-meson couplings. 



IV. FINITE TEMPERATURE RESULTS 



Taking the parameters in Eq. (|3ip . one can solve Eqs. Q to calculate the mean field values 
a u , a s and 4>3 at finite temperature. The behavior of effective quark masses and condensates, as 
well as the curves for the traced Polyakov loop <£, are similar to those obtained in Ref. [lgI within 
an ILM-motivated nonlocal chiral model. The discussion of those results is qualitatively the same 
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as in our case, therefore it will not be repeated here. We just state that, as expected, there is 
a crossover phase transition in which chiral symmetry is restored, and consequently one finds a 
sharp peak in the chiral susceptibility. The transition temperature (defined as the position of this 
peak) is found to be T c = 202 MeV. This value is in much better agreement with lattice results, 



r patt) = 16Q _ 20Q MeV ^ than the value recently oi3 tairiec i in tne i oca i gu(3) PNJL model, 



jt(pnjl) _ M e y j3g[ | j n addition one finds a deconfinement phase transition, which occurs at 
about the same critical temperature. 

We concentrate here in the evolution of meson masses and decay constants with temperature, 
which has not been previously addressed in the context of SU(3) nonlocal models. Pseudoscalar 
meson masses can be determined by solving Eqs. (|2ip . while the same procedure applies to the 
scalar meson sector replacing upper indices "— by "+" in Eqs. (|18H19p . As discussed in Sect. II, 
these values correspond to the spatial screening-masses of the mesons' zeroth Matsubara modes. 
Our numerical results are shown in Fig. 1, where we quote the values of meson masses as functions 
of the temperature. In Fig. 1(a) we show the behavior of the pseudoscalar mesons it and r] together 
with the curves for the scalar mesons a and ao, which are chiral partners of the former. It is seen 
that pseudoscalar meson masses remain approximately constant up to the critical temperature 
(this is reasonable, since they are protected from chiral symmetry), while scalar meson masses 
begin to drop at about 150 MeV. Above T c pseudoscalar masses get increased, in such a way 
that they become degenerate with the masses of their chiral partners, as expected from chiral 
restoration. In particular, the fact that this occurs right after the transition in the case of the 
(77, ao) pair indicates that the strange contents of the rj and ao mesons become suppressed above 
the critical temperature. When the temperature is further increased, all four masses are found 
to rise continuously, showing that now the mass is basically dominated by thermal energy. At 
very large temperatures the curves should approach asymptotically the value corresponding to a 



qq pair of uncorrelated massless quarks rn^ = 2nT 36(. At T ~ 300 MeV, however, the Polyakov 
loop has not yet reached its asymptotic value 4>^/T\t^oo = 0, and it still provides a non-negligible 
correction to the quark screening mass. In fact, we find ^3/^1^=300 — 0.93. Thus, around this 
temperature we expect rn^j = 2(irT — fo), which is shown by the dotted lines in Figs. 1(a), (b) and 
(c). For vanishing quark dynamical masses, this value of tum corresponds to a pole in the n = 



mode for the integrals in the functions C^(— k 2 , 0), see Eq. flT3J). Indeed, as discussed in Ref. 



Eqs. (I21|) can be satisfied only in the vicinities of these poles. On the other hand, in general it is 
seen that the functions C^(—k 2 ,0) [and therefore also the functions Gm(— £; 2 ,0)] are well defined 
for low values of k. If k is increased, at some critical point k CI it, usually called "pinch point", 
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1500 





FIG. 1: Solid lines quote scalar and pseudoscalar meson masses as functions of the temperature. Dotted 
and dashed-dotted lines stand for the value = 2(irT — fa) and the qq production threshold, respectively 
(see text). 
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the integrals become divergent and need some regularization 



prescription. In the present work we 



follow the prescription discussed in the Appendix of Ref. 23] , conveniently extended to the finite 
temperature case. The pinch point occurs when both effective quarks are simultaneously on-shell, 
thus it can be interpreted as a threshold above which mesons could decay into two massive quarks. 
In Fig. 1(a) this threshold is represented with the dashed-dotted curve (above T c , it approximately 
matches the value of rn3 mentioned previously). It can be seen that all four meson masses in 
Fig. 1(a) remain below the threshold for the temperature range considered. 

In Fig. 1(b) we represent the curves for the masses of the pseudoscalar mesons K, and their scalar 
partners k. It is seen that for some temperature range the equation G K {— k 2 , 0) = has no solution 
for real k, therefore the mass is defined as the minimum of the function G K (—k 2 ,0), as discussed 
in the previous section. These mass values correspond to the dashed stretch of the corresponding 
curve. It is worth to notice that the K and k meson masses match only at T ~ 225 MeV, i.e. at 
a temperature somewhat larger than T c . This is clearly a consequence of the large current strange 
quark mass, which is expected to move the SU(3) chiral restoration to higher temperatures. Finally, 
in Fig. 1(c) we quote the temperature dependence of fo and rf masses. As before, dashed stretches 
in the curves indicate the regions in which the corresponding function GM(—k 2 ,0) has no zero 
for real k and, therefore, the mass uim is defined by the position of its minimum. Let us first 
focus on the behavior of the rj mass. In contrast with some results found in Ref. 38j, where the 
corresponding temperature dependence has been studied in the framework of a Dyson-Schwinger 
approach, we do not observe any kind of enhancement of m„/ around T c . It should be noticed 



that in the framework of Ref. 



38 ] the effect of the U(l)yi anomaly is modelled in a simpler way, 



namely by considering it only at the level of mass shifts. In this sense our result is consistent 
with the analyses of the rf pole mass performed within the local SU(3) PNJL model [35] and the 



quark-meson model 



39], where no enhancement was found either. Concerning the degeneracy of 



rf with its chiral partner fo, we see that such a degeneracy is achieved only at T ~ 300 MeV. This 
a consequence of the strange quark contents of these mesons, which, as we will see below, become 
larger as the temperature increases. 

Next, in Fig. 2 we quote the behavior of pseudoscalar meson decay constants, which can be 
calculated from Eqs. (f26j) . ([26]) and (j28|) . Fig. 2(a) shows the curves corresponding to the decay 
constants f n and fx- It is seen that both decay constants drop at the phase transition. We observe, 
however, that due to the strange quark content of the kaon the corresponding decay constant shows 
a slower decrease after the transition. The behavior of the decay constants associated with r} — rf 
system is shown in Fig. 2(b). In the case of f VQ and f Vs such a behavior is similar to that of f n , 
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FIG. 2: Behavior of pseudoscalar meson decay constants as functions of the temperature. 




FIG. 3: (a) Behavior of meson mixing angles as functions of the temperature. The dashed line shows 
the "ideal" mixing angle #idcai = tan -1 \/2. (b) Behavior of the topological susceptibility relative to its 
T = value as function of the temperature in the nonlocal model (solid) and in the PNJL SU(3) model 
3 (dotted). 

while the decrease after the transition is less pronounced for f^i . Again, this behaviour of the rj 
decay constants can be understood in terms of its larger strange quark content. Here we have left 
blank the range in which the rf mass is not well denned. 

In Fig. 3(a) we plot the behavior of the mixing angles 6^ and 9„>, which can be calculated from 
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FIG. 4: Behavior of quark-meson couplings as functions of the temperature. 

Eq. (I22D . It is seen that above the phase transition both angles tend to a common value, which is 
natural since meson masses also tend to unify. More interestingly, they converge to the so-called 
"ideal" mixing angle #;deai = tan -1 y/2 ~ 54.7° (dashed line in the figure). This means that, as 
suggested above, the rj meson becomes approximately non-strange, while rj approaches to an ss 
pair. The same happens with the a — fo pair (in the figure we have quoted only 9 a , since the 
fo meson mass lies above the qq threshold). The fact that the mixing angles go to the "ideal" 
value for large temperatures implies that the U(1)a anomaly tends to vanish in this limit. Another 
signature of this fact is that axial chiral partners (ir, rf) and (cr, ao) become almost degenerate at 
T ~ 300 MeV. However, perhaps the best indication of the vanishing of the U(l)^ anomaly is 
provided by the topological susceptibility x which, in pure color SU(3) theory, is related to the r/, 
r] and K masses through the Witten-Veneziano formula 

6 



f2 ^ 



2 , 2 



2m\ 



(33) 



Various existing lattice calculations 



401 ] show a sharp decrease of \ a t the critical temperature. In 



our framework, the topological susceptibility can be calculated from 



H 2 



X 



2C UU S U S S -\-C ss S u 



S U {S U + 2S S ) i <■ iks 
Su(S u + 2S S 



V2 




(34) 



y 00 / \ ^80 

where Q is a 2 x 2 matrix whose matrix elements are given in Eq. (I19p . and all functions are evaluated 
at (q 2 ,^) = (0,0). This expression has been obtained following similar steps as those described 
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in Ref. |4l| for the case of the (local) NJL model. The numerical results are shown in Fig. 3(b). To 
be able to compare with the result obtained in the local PNJL SU(3) model 35j (where, as already 
mentioned, T c turns out to be too high), we show the normalized value of x 1 ^ 4 as a function of 
T/T c . For both models one finds a sharp decrease in the topological susceptibility at the critical 
temperature, this decrease being steeper in the nonlocal model. Indeed, at T/T c = 1.5 the ratio 



X 



1/4/ 



X 



1/4 
(T=0) 



is about 11% for the local PNJL model, while for the nonlocal model it is roughly 

1/4 



one half of this value. The value of Xrp—o) * s f° un d to be about 162 MeV in the nonlocal model, 
while one gets ~ 180 MeV in the PNJL. Recent lattice calculations (see Ref. 42J and references 



therein) indicate that x\t=o) ~ ^0 MeV in pure gauge theories. However, light dynamical quarks 

43J. For example, in the lattice calculation 



are expected to suppress the topological susceptibility 



carried out in Ref. 



44 ] the authors find x 



1/4 
(T=0) 



163 MeV for a two-flavor case in the region where 



the current quark masses are around 20 MeV. 

For completeness, we conclude our description by quoting in Fig. 4 the behavior of the quark- 
meson couplings GMq- These can be calculated from Eqs. (|23p and (|24p for pseudoscalar mesons, 
and similar relations hold for the scalar meson sector with the appropriate changes in the functions 

cv,0f 2 ,o). 



V. SUMMARY AND CONCLUSIONS 



In the present work we have studied the finite temperature behavior of light scalar and pseu- 
doscalar meson properties in the context of three-flavor nonlocal chiral models that include mixing 
with active strangeness degrees of freedom. The effect of gauge interactions has been introduced 
by coupling the quarks with a background gauge field, and the deconfinement transition has been 
studied through the behavior of the traced Polyakov loop. For a given parameterization of the non- 
locality — which, for simplicity, here is introduced through an exponential form factor — , at zero 
temperature the model has five free parameters. We have chosen to fix the average non-strange 
quark mass m to a phenomenologically sound value of fh = 5 MeV, whereas the remaining four 
parameters have been determined by requiring that the model reproduces correctly the measured 
values of the masses of the pion, kaon and rf pseudoscalar mesons, and the pion decay constant 
f n . Using this set of parameters one can obtain a very good description of the remaining zero tem- 
perature pseudoscalar meson properties, as well as adequate values for the scalar meson masses. 

In the extension to finite temperature the former parameter values have been kept fixed, while 
those appearing in the Polyakov loop potential have been taken from a fit to lattice results. As 
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expected, the model shows a fast crossover phase transition, corresponding to the restoration of 
SU(2) chiral symmetry. The transition temperature (defined as the position of the peak of the 
corresponding chiral susceptibility) is found to be T c = 202 MeV. This value is in better agreement 
with lattice results, namely J^ latt ^ = 160 — 200 MeV [3J], than the value recently obtained in the 



local SU(3) PNJL model, T C (PNJL) = 259 MeV 



351 ]. In addition one finds a deconfinement phase 



transition, which occurs at about the same critical temperature. Concerning the behavior of meson 
masses with temperature, it is seen that pseudoscalar meson masses remain approximately constant 
up to T c , while scalar meson masses begin to drop at about 150 MeV. Beyond T c pseudoscalar 
masses get increased, in such a way that they become degenerate with the masses of their chiral 
partners, as expected from chiral restoration. The temperature at which chiral partners meet 
depend on the strange quark composition of the corresponding mesons, i.e. the masses of mesons 
containing no strange quarks match almost immediately after T c , while fo and rj masses meet only 
at about 1.5 T c , the situation being intermediate for K and k mesons. Regarding the properties 
of the rj — rj' sector, it is seen that the corresponding mixing angles tend to converge to the so- 
called "ideal" mixing, which indicates that the 11(1)^ anomaly tends to vanish as the temperature 
increases. This is also seen in the behavior of the topological susceptibility which, as expected 
from lattice calculations, shows a sharp decrease after the chiral phase transition. It should be 
noticed, however, that in the present nonlocal model such a decrease is faster than that obtained 
in the local PNJL SU(3) [35]. Finally, we notice that, in agreement with the local model — and in 
contrast with what was suggested in the framework of a Dyson-Schwinger approach [38 ] — we do 
not observe any kind of enhancement of the rj mass around the critical temperature. 
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